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We present a study of two-photon decays of light mesons, focusing on the neutral pion decay. This 
important process highlights the effects of the axial anomaly in QCD but has been little studied 
on the lattice. By applying the Lehmann-Symanzik-Zimmermann (LSZ) reduction formula, we 
reconstruct the electromagnetic matrix elements from three-point vector-vector Green functions 
calculated on 2h-1 -flavor isotropic clover lattices. 
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1. Introduction 

The decay of the neutral pion into two photons has attracted a great deal of scrutiny since the 
early days of particle theory. In the Standard Model, this decay is understood to be primarily due 
to the effects of the Adler-Bell-Jakiw axial anomaly, where it is much larger than any competing 
process. Despite its long history, the 71° — > 77 decay is difficult to study experimentally, since both 
the initial and final states are electrically neutral. The large uncertainty in this width (7.84(56) eV, 
about 7%) propagates into determinations of the splitting between the up- and down-quark masses. 

Given its importance as a fundamental prediction of the Standard Model and the relatively poor 
precision with which it is currently known, it is not too surprising that there is a substantial experi- 
mental program attempting to improve this measurement. The Primakoff Experiment (PrimEx) in 
Hall B at Jefferson Lab (JLab) uses photon fusion between a real photon and the electric field of a 
heavy nucleus to produce neutral pions. Combined with a high-precision calorimeter, PrimEx was 
constructed to provide a percent-level precision measurement of the neutral pion width. 

From the theoretical side, the situation is somewhat better. The amplitude of the neutral pion to 
two photons process is exactly given by the axial anomaly in the chiral limit, and chiral perturbation 
theory (XPT) may be used to give a physical result. To leading order, XPT gives 

which is in apparent agreement with experiment. However, there are some complications. At 
next-to-leading order, XPT must consider the mixing of the pion with the eta and eta-prime. This 
increases the prediction to 8.10(8) eV[l]. Lattice QCD should also be able to determine this quan- 
tity from first principles, but since XPT appears to perform so well, this will simply be a crosscheck. 
On the other hand, lattice calculations can address other related quantities where XPT is expected 
to perform poorly, such as the processes T] ^ 77, T]' —> 77 and ao — > 77- Measurements of these 
quantities will be performed by PrimEx following the 12 GeV upgrade at JLab. In addition, ex- 
periments using photoproduction (such as GlueX) will probe photon-fusion production of neutral 
mesons (77* <I>°). 

Despite the clear motivations for making a measurement of two-photon decays of neutral 
mesons, lattice calculations have done little to address the issue. The difficulty faced in lattice 
QCD is simple: the photon is not an eigenstate of QCD. An operator carrying the quantum num- 
bers of the photon in lattice QCD will instead create a rho meson (or two pions, depending on the 
quark mass used). The resolution of this difficulty, proposed by Ji and Jung[2, 3], uses perturbative 
field theory techniques to express the photon as a superposition of QCD eigenstates accessible on 
the lattice. This technique was recently used to study the two-photon decays of chaimonium in 
Ref. [4]. 

2. Lehmann-Symanzik-Zimmermann Reduction 

The method of Ji and Jung applies the Lehmann-Symanzik-Zimmermann (LSZ) reduction to 
evaluate a photonic matrix element. We are interested in matrix elements of the form 

{y{qiMy{q2,h)\^{p)). (2.1) 
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where the two photons have momenta ^12 and polarizations Ai,2 and the neutral meson 4> has 
momentum p. Using a standard LSZ reduction, we re-express this using a time-ordered (for now, 
working in Minkowski spacetime) product of photon fields: 

- lim el!>iquXi)ei^>iq2,hWi^l2 ^d'xd'ye''^^■y+'^2■^o\T{A^'iy)A^x)}Mp)), (2.2) 

where e is a polarization tensor and A is the vector photon field, which we have Fourier transformed. 

Since lattice QCD can only directly measure eigenstates of QCD, we cannot treat these A 
fields. One solution to this might be to include explicit QED U(l) fields during our lattice gauge 
generation, but since the cost of creating new ensembles is so great, we would rather attempt to 
express them in terms of quark and gluon operators accessible to QCD alone. To do this, we apply 
perturbative QED to integrate out the photon fields: 



DAD\ffD\i/e'^'i''''A^{y)A''{x) «i 

jDADYDxife''' (^... + ^ [i^/fiMp] {z)[vfvA^]{w) + ..)jA^{y)A''{x), (2.3) 

where we take a path integral over all possible configurations of the photon field A and the quark 
fields Y and iff (with electric charge e) using the free-field action 5o and the interacting QED action 
Sqed- This can be simplified by applying Wick contraction to the photon fields: 



lim el!>ei^>q'M' ^d'xd'wdW^^■''D^'P{0,z)D''^{x,w){0\T{jp{z)ja{w)}Mp)), (2.4) 

where j^l = xj/Yn Xj/ and D is the photon propagator, 

«'"(0,.) = -,.'"/(4^. (2.5, 

By substituting the expUcit form of the photon propagator, we cancel the inverse propagators 
and most of the integrals over spacetime become simple delta-functions. The remaining expression 
is 

ehl'^*e^^^* Jd'ye'''-y{0\T{f{0)f{y)}Mp)). (2.6) 

We now need to rotate into Euclidean space. This is possible as long as our integration contour 
does not hit a pole, where one of the photons can mix with an on-shell particle. Therefore, we must 
keep < Mp (or E'^„, depending on the quark mass). 

The resulting expression is something that can be evaluated in lattice QCD: 



2E^(p) 



yjdtie-'^^^''-'^ (^r Ud'ze-'P-^'Mltf) Jd'ye'^^'yriy,t)fid,ti)\\, (2.7) 



where Wi is the energy of the first photon, (p^ is a creation operator for the neutral meson and 
and are its overlap factor and energy respectively. The expression between the angled brackets 
is just the three-point correlation function with a meson on one end and vector currents at the other 
end and inserted. The remaining parts describe how to combine these QCD currents into a photon 
of the appropriate energy. 
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3. Lattice Calculation 



The most straightforward way to evaluate Eq. 2.7 on the lattice is to compute the three-point 
function with all insertion times and perform the integral by explicit summation. For this first test 
of the light-meson two-photon decay calculation, we use lattices made available through the ILDG 
by CP-PACS. Particularly, we use 20^ x 40 lattices with 2-1-1 dynamical flavors of clover fermions 
having pion mass Mj^ = 725 MeV and lattice scale c?^' 2.25 GeV[5]. 

For the explicit-integral ("slow") method, we fix the location of the neutral meson at timeslice 
tf = 32 and zero momentum. We generate propagators using the same clover action used in the 
dynamical sector under Dirichlet boundary conditions and apply Gaussian smearing to improve 
overlap with the ground-state mesons. Propagators ai^e generated from point sources at tf, from 
these, we create sequential- source propagators applying the smearing and momentum-projection 
at tf. These propagators are then combined into three-point functions by tying them together with 
various momentum projections (0 < |^2p < 5) at timeslice t and including the appropriate gamma 
structure for vector currents. (See Figure. 1 .) 





Figure 1: Left: Schematic of the "slow" method. The black arrow indicates an ordinary propagator; the 
blue arrow is a sequential-source propagator. The three-point function must be evaluated for all f, and t. ^ 
denotes momentum projection. Right: Schematic of the "fast" method. The sequential propagator source 
has support on all timeslices and is weighted by the exponential factor. Note that the propagator source is 
now at tf, reversed relative to the slow method. 



In order to ascertain whether our calculation will suffer from lattice distortions, we need to 
scrutinize the time-dependence of the integrand in Eq. 2.7, 



(3.1) 



where '^pvv is the three-point correlator. If the integrand is peaked too sharply, we will not be 
able integrate it accurately; if it is too wide, we cannot capture the integral within the lattice time 
extent. In addition, we need to check for distortion due to the boundaries and proximity to the sink 
timeslice. A sample of our results is shown on the left-hand side of Figure 2. The peak is well 
resolved, neither too nanow nor too wide; it becomes cut off by the ?, = boundary for the red 
f = 8 peak and is seriously distorted by the tf = 32 sink for the brown t = 32 peak. We must avoid 
these extremes when looking for a plateau in the integrated matrix element. 

We perform the integral explicitly by summing over all f, for each value of t. The integral will 
be nonzero only when e^vpa£^£^(]^i<l2 / 0. A sample plateau is shown on the right-hand side of 
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Figure 2: Left: Dependence of the integrand on selected t, denoted by color, and f,, along the x-axis. 
Right: Integrated matrix element for x, {/i, v} = {^jZ}- Fitted value over the plateau region shown as a 
pink bar 



Figure 2. As previously discussed, we expect to see a plateau away from the boundary at ? = 
and the sink at ? = tf; we choose to fit to the range t G [10,25]. Although it is possible for there to 
be exponential contamination from excited states, no such problem is seen here. This may be due 
to the large gap between the ground-state pion and its first excited state or due to the excited pion 
having small coupling to the two-photon state. 

We calculate the neutral pion decay for a variety of photon virtualities. Due to the slow-method 
construction we use, the value of (0\ (or equivalently —Q\ = (of — \q\ may be set arbitrarily. The 
value of Q\ is then determined by conservation of momentum and energy. We plot the decay 
constant as a function of Q\ and Q\ on the left-hand side of Figure 3. The data are well described 
by a monopole fit: 

where F and Mpoie ai^e fit independently for each Q\. Vector-meson dominance (VMD) suggests 
that the pole mass should be approximately equal to the rho mass. We plot M^^^^ as a function of 
Q\ on the right-hand side of Figure 3. The agreement with the lattice rho-meson mass is relatively 
good. 

If we assume VMD and set the pole mass to Mp , we can use a simpler form: 

where F is allowed to vary with Q\. We restrict the fit to regions well away from the pole to 
avoid distortion. This fit is shown on the left-hand side of Figure 4. 

Since the two photons are indistinguishable, there should be symmetry under exchange of Qy 
and Q2. We therefore try fitting all the data simultaneously to a double-pole form: 

/ + e?/M2„,e)(l+el/^pole 

where now there is only one F and Mpoie. This fit is shown on the right-hand side of Figure 4. 
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Figure 3: Left: The neutral meson two-photon decay matrix element as a function of and Q\. Different 
values of Q\ (-0.5 GeV^ to 1.0 GeV^ in steps of 0.25 GeV^) are displayed as different colors. ^ runs along 
the X-axis. The bands show the fit according to the form in Eq. 3.2. Right: Pole mass as a function of Q\, 
extracted from fit 




GzCGeV^) 



Figure 4: As Figure 3 with fits according to Left: Eq. 3.3 and Right: Eq. 3.4. Only data included in the fit 
are shown. 



This small data sample shows that our method for extracting the matrix element of the neutral 
meson decay to two photons is working. However, the technique is very expensive, since we 
must calculate T (the time-extent of the lattice) different propagators to sample the integrand. 
We can greatly improve this by using a sequential source trick, shown on the right-hand side of 
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Figure 1. Rather than calculating propagators with different ?, individually, we take the solution to 
the first propagator, multiply it by the exponential factor g^'^Kf--') and use this whole quantity as a 
sequential source. The resulting sequential propagator implicitly contains the integration. 

This "fast" method has a few disadvantages: We can no longer directly examine the integrand 
to be sure that we are con^ectly evaluating the integral. Since we have already invested in a test set 
using the slow method, this is not a problem. We can also no longer vary Q\ without recalculating 
the sequential propagator, since C0\ appears in its source. However, since we save a factor of T , 
which can be very large, the fast method will be our method of choice in future work. 

4. Summary 

We have demonstrated that the method of Jung and Ji allows us to access the two-photon 
decays of neutral mesons in lattice QCD. Using this method on two-flavor lattices, we confirm the 
predictions of vector-meson dominance, that the matrix element has a monopole form with pole 
mass approximately equal to the rho mass. In the future, we will apply the fast method to greatly 
expand the range of quantities we calculate, including conserved-current insertions and scalar and 
axial meson decays. 
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